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Abstract After a brief recapitulation of the general interest of parton densities, we discuss multiple 
hard interactions and multiparton distributions. We report on recent theoretical progress in their QCD 
description, on outstanding conceptual problems and on possibilities to use multiparton distributions 
as a laboratory to test and improve our understanding of hadron structure. 
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1 Parton distributions: a brief recapitulation 

Among the most interesting and most difficult aspects of QCD is the relationship between quarks and 
gluons, which are the basic degrees of freedom of the theory, and hadrons, which are the physical states 
observed at macroscopic distances. Parton distributions are among the most prominent quantities that 
describe this relationship. Important issues in this context are: 

— Parton densities describe quarks, antiquarks and gluons as they manifest themselves in short- 
distance processes, and they have a direct connection to the fields that appear in the QCD La- 
grangian. How are these degrees of freedom related with the quarks that appear in non-perturbative 
approaches to hadron structure, such as constituent quark models, chiral quark models, the Dyson- 
Schwinger approach etc.? 

— The appearance of partons at small momentum fraction x can to a large extent be described as the 
result of parton radiation in the perturbative regime, as it is for instance encoded in the DGLAP 
equations. However, in particular the fits performed by the Dortmund group clearly show that 
perturbative physics alone is not sufficient for understanding sea quarks and gluons, given that they 
must already be present at low resolution scales, where they can only be of non-perturbative origin. 
Distributions such as u — d ot s — s are particularly clear indicators of this fact. 

— Which roles do confinement and chiral symmetry breaking play in determining the distribution of 
partons inside a hadron? 

— Gluons and the choice of gauge play an essential role when one describes the dynamics of partons 
inside hadrons. Our understanding of this issue has seen important progress in the last decade but 
is still far from satisfactory, as we will see at the end of this section. 

Parton distributions come in different varieties, which address different aspects of hadron structure. 
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PDFs, i.e. parton distribution functions /(a:), are the familiar functions describing the distribution of 
partons in longitudinal momentum fraction x. There is a vast number of processes where they can 
be measured. The PDFs for unpolarized partons in the proton are an indispensable ingredient for 
understanding high-energy lepton-hadron and hadron-hadron collisions, and with the exception of 
certain distributions they are among the most precisely measured non-perturbative quantities in 
QCD. 

TMDs, i.e. transverse-momentum dependent distributions /(x,k), describe the joint distribution of 
partons in their longitudinal momentum fraction and their transverse momentum. They are relevant 
in processes with a measured transverse momentum much smaller than the hard scale, e.g. in Drell- 
Yan production when the transverse momentum of the lepton pair is much smaller than its invariant 
mass. (Here and in the following, vectors in the transverse plane are written in boldface.) 
TMDs quantify a variety of spin-orbit correlations at the parton level. A particularly intriguing 
aspect of their dynamics is the role of gluons that are described by Wilson lines (see below). We 
note that the term TMD also includes transverse-momentum dependent fragmentation functions. 

GPDs, i.e. generalized parton distributions F{x,^,t), appear in the description of exclusive processes 
like deeply virtual Compton scattering or meson production. The Mandelstam variable t can be 
traded for the transverse momentum difference A between the incident and the scattered proton, 
which is Fourier conjugate to the transverse position b of the parton inside the proton. If the 
skewness variable ^ (which describes the longitudinal momentum transfer to the target) is zero, the 
Fourier transformed distributions f{x, b) describe the joint density of partons in their longitudinal 
momentum fraction and their transverse position 3] . One often refers to /(a;, b) as impact parameter 
distributions. 

The Mellin moments J dx x"'^^ F{x, ^,t) of GPDs can be identified with the form factors of local 
operators. This provides in particular a connection between GPDs and the electromagnetic Dirac 
and Pauli form factors of the nucleon, which have been measured with great precision over a wide 
range of t. The moments of GPDs can also be evaluated in lattice QCD [4J. 

A common feature of these different types of distributions is that they are defined in terms of matrix 
elements 

{P2\qi3{z2)W{z2,Zi)qa{zi)\pi) (1) 

of quark-antiquark operators between proton states. For gluon distributions, the fields q and q are 
replaced by the gluon field strength These matrix elements have a rather simple interpretation 
in light-cone quantization, where the fields 9(22) and q{zi) at zero light-cone time can be expanded in 
terms of creation and annihilation operators for quarks or antiquarks (see e.g. (H). A subtle point of 
interpretation arises due to the Wilson line W{z2, zi), which will be discussed shortly. 

Longitudinal and transverse coordinates play quite different roles in the matrix element ([1]), which 
means that one has lost manifest three-dimensional rotation invariance even if both proton states pi and 
P2 are at rest. This corresponds to the fact that in processes where parton distributions can be observed, 
there is a physically preferred direction. The variables of and in ([1]) are related to the light-cone 
plus-momenta of the quark or antiquark by a Fourier transform, and a parton interpretation is most 
straightforward in a reference frame where the protons move fast, i.e. where and are large. The 
transverse variables in ((1} differ between PDFs, TMDs and GPDs. One can switch between transverse 
position and transverse momentum of the proton states or the partons by a Fourier transform. There is 
a subtle point in this, which can easily be explained by a simple calculation. The Fourier transformed 
quark field 



g(z-,k) = J d^z e"'^q{z-,z) (2) 

at z^ = annihilates quarks with transverse momentum k and creates antiquarks with transverse 
momentum — k. For a bilinear operator as in ([T]) we then have 

g-(k) g(k') = J d'z d\' e-(k-k'z') ^ 
where here and below we suppress plus- and minus-coordinates for simplicity. With 

kz - k'z' = (z - z') + (k - k') ^ (4) 



a 
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Fig. 1 Graphs with gluon exchange between hadrons and the hard-scattering subprocess for semi-inclusive 
deep inelastic scattering (a), Drell-Yan production (b) and hadron pair production in pp collisions (c). 



we see that in the region where q creates and q annihilates a quark, the average transverse momentum of 
the created and annihilated quark is Fourier conjugate to the difference of their transverse positions, 
whereas the difference of the transverse momenta is conjugate to their average transverse position. 
Taking the matrix element of ^ between proton states that are localized at transverse position 
(z -I- z')/2 — b, and making an appropriate Fourier transformed with respect to the minus-positions 
of the fields, one obtains Wigner distributions M^(k, b) that depend on both the average transverse 
momentum of a quark and its average transverse position relative to the proton, where the "average" 
refers to the two fields in the operator ([3]). The integral / cPh VF(k, b) gives the density of quarks with 
transverse momentum k, i.e. a TMD, whereas the integral / (i^kW^(k, b) gives the density of quarks 
with transverse position b in the proton, i.e. essentially a Fourier transformed GPD. This shows 
explicitly that TMDs and GPDs contain complementary information about partons in the proton, and 
that both types of quantities are descendants of higher-level functions, namely of Wigner distributions. 
For a general introduction to Wigner distributions we refer to . 

Let us now return to the Wilson line in (jlj, which we have glossed over so far. When establishing 
factorization for processes such as Drell-Yan production or semi-inclusive deep inelastic scattering 
(SIDIS), one finds that an arbitrary number of gluons in a right-moving hadron can attach to the 
hard-scattering process without any twist suppression, provided that the gluons have polarization in 
the plus-direction. This is illustrated in Fig. [T^ and b. It is the effect of these gluons that is summed 
up in the Wilson line W{z2, zi), with a path between zi and Z2 that depends on the process. 

For k integrated distributions, the situation simplifies if one works in the light-cone gauge = 0, 
where the Wilson line reduces to unity so that one has a literal parton model interpretation of the 
operator in (jlj as creating and annihilating quarks or antiquarks. For k dependent distributions, 
however, there are subtle spin effects such as the Sivers or Boer-Mulders asymmetries, which are quite 
explicit in a covariant gauge, whereas in = gauge they are due to Wilson line pieces at infinity, 
corresponding to gluon modes with /c+ = 0, see e.g. [3;|lj. 

Closer inspection shows that in TMDs one can actually not take the Wilson line W{z2, zi) along a 
path in the light-like direction, since this leads to rapidity divergences @. To regulate these divergences, 
one can introduce a direction v with nonzero plus- and minus-components, so that the resulting Wilson 
lines reduce to unity in the axial gauge Av = rather than in light-cone gauge. A formulation of 
factorization with TMDs in the framework of light-cone quantization and the — gauge has so far 
not been achieved. On the positive side, the formulation with Wilson lines along a direction v allows 
one to resum Sudakov logarithms to all orders. This has led to a successful phenomenology of processes 
with small measured transverse momentum, see [l^ and references therein. On a more fundamental 
level, the physical phenomena associated with TMDs force (and guide) us to think more deeply about 
the role of gluons and of gauge degrees of freedom when describing how parton emerge from hadrons 
and interact with each other. 

Finally, we must note that TMD factorization has so far only been established for processes with a 
simple color structure in the hard scattering, such as SIDIS and Drell-Yan production. Serious obstacles 
for such a formulation have been identified for processes like the production of back-to-back hadrons 
or jets in pp collisions, where graphs like the one in Fig.[T]: do not appear to admit a resummation into 
Wilson line operators It remains an open question whether a factorized description (and along 
with it, Sudakov resummation with the same accuracy as for Drell-Yan production) can be found for 
this important class of processes. 



Fig. 2 Double hard scattering processes producing a Drell-Yan lepton pair and a dijet (a) or two Drell-Yan 
pairs (b). 

2 Multiparton interactions: what are they and why are they interesting? 

In a generic hadron-hadron collision, several partons in one hadron scatter on corresponding partons 
in the other hadron. At high energy, several of these scatters can have a hard scale and produce 
particles with large invariant mass or large transverse momenta. Examples are shown in Fig. [31 The 
effects of such multiparton interactions (also referred to as multiple interactions) average out or are 
power suppressed in sufficiently inclusive observables. This is why they are not included in the familiar 
factorization formulae for hadron-hadron collisions, which involve single-parton densities and a single 
hard-scattering subprocess. However, for more exclusive observables and more specific kinematics, 
contributions from multiple interactions can be substantial, and it has been estimated that they will 
play an important role in many analyses at the LHC. There is a long history of theoretical and 
experimental studies of such interactions, and the field has seen a strong boost in activity since the 
LHC has started operation, see e.g. the recent proceedings [T2l |. 

The description of multiparton interactions involves multiparton distributions, which are interesting 
in their own right since they contain information about correlations between different partons in the 
proton wave function. This information is inaccessible in the singe-parton distributions we discussed 
in Sect. □ 

The existing phenomenology of multiple interactions (including implementations in Monte Carlo 
event generators) is based on a simple and physically intuitive picture, which involves however many 
simplifying assumptions. A systematic treatment in QCD, which would match the level of sophistica- 
tion achieved for single hard-scattering processes, remains to be given. Some steps into this direction 
have been taken in |l3|. and in the present proceedings we will point out results that have been ob- 
tained, open problems that have been identified, and prospects of using multiple interactions to learn 
more about hadron structure at the parton level and about theoretical approaches to understand this 
structure. 



3 Some basic results 

To begin with, let us specify the theoretical framework we will be using. We consider double hard 
scattering as the simplest (and often most prominent) case of multiparton interactions. Requiring both 
parton-level scattering processes to be hard allows us to use the concept of factorization and the predic- 
tive power of perturbation theory. Since the phenomenological interest in multiple interactions comes 
from the necessity to understand the hadronic final state in detail, we keep the transverse momenta of 
the produced particles differential and are hence led to using transverse-momentum dependent factor- 
ization and the associated multiparton distributions. Given the status of TMD factorization for single 
hard scattering described above, we consider the double Drell-Yan process of Fig. as a case where 
the prospects for obtaining rigorous theory results are best, referring the practically important cases 
where secondary interactions produce jets to future work. 

The graph in Fig. [3^ shows the production of a gauge boson pair by two hard-scattering processes. 
This graph can be evaluated with the same techniques that are used for single hard scattering, where 
small momentum components are neglected compared to large ones in each parton-level subprocess. In 
this approximation, the plus-momenta of right-moving partons and the minus-momenta of left-moving 
ones are fixed by the final-state kinematics, but their transverse momenta are not. With the momentum 
assignments in Fig. one finds = qf /p^ and = q~ /p~ for the longitudinal parton momentum 



a 
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Fig. 3 Graphs for the production of two electroweak gauge bosons by double (a) and by single (b) hard 
scattering. The dashed line indicates the final-state cut, and the decay of the bosons into leptons is not shown 
for simplicity. 



fractions. Momentum conservation gives r + f = for the transverse components of the momentum 
mismatch between partons on the left and the right of the final-state cut. The contribution of Fig.|3^ to 
the cross section is then proportional to / cPr Fai,a2 i^ii kj, r) Fa^ a^ {xi, k^, — r), where the first factor is 
the distribution for two quarks and the second factor the distribution for two antiquarks in the proton. 
Fourier transforming the distributions w.r.t. their last argument, one finds that the conjugate position 
variables satisfy y = y, and the full expression of the cross section reads 



da 



Fig. [Hi 



Yii^i dxi dxi d'^qi 



X (Ti,aiai(9i) 'T2,a2a2(92) / c'^y ^ai ,03 (a;^ , , y ) (xi , k, , y ) 



(5) 

where ^i^aiai denotes the hard-scattering cross section for single-boson production. The statistical 
factor 5 is 2 if the produced bosons are identical and 1 if they are not. 

From our discussion in Sect. [I] it follows that the Fourier conjugate variable y of r represents 
the distance between the two scattering partons, averaged between the scattering amplitude and its 
complex conjugate. On the other hand, ki and k2 are the transverse momenta of the partons, averaged 
in the same sense. Faj^a2 (2:4, k^, y) thus has the structure of a Wigner distribution in the transverse 
degrees of freedom. It is gratifying that this allows for a rather intuitive interpretation of the cross 
section formula we have two quark- antiquark annihilation processes separated by a transverse 
distance y, where in each annihilation the incident transverse parton momenta k^ and k^ add up 
to the measured transverse momentum of the produced gauge boson. Note, however, that the 
result ^ is fully quantum mechanical and has been obtained from Feynman graphs using standard 
approximations, without the need to appeal to semi-classical arguments. The operator definition of a 
two-quark distribution reads 



-^01,02 (•^i: ki, y) 



n 



2p+ dy-{p I g(-iz2) Wi-^Z2, \z2) 9(5^2) 



X 9(2/ - ^{V -\z\,y+ Fa, q{y + \zi) I p) 



(6) 



and is a natural extension of the definition of TMDs for a single parton described in Section [TJ 

The transverse- momentum integrated distribution ,02 (a^i I y) = / rf^ki (i^k2 Faj^a2(xi, k^, y) gives 
the joint density of two quarks with momentum fractions X\ and xi and relative transverse distance y. 
It enters in the transverse-momentum integrated cross section as 



da 



Fig. Hi 
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O'l ,a2— QiZiiji^g 
O'l ,a2— g,/ig,<5g 



O-l.aiSilg?) CT2,Q2a2(92) J d'^yFa,, 



,Xx^,y) Fai,a2(a;j,y) 



(7) 
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The integration over ki and k2 puts the relative transverse coordinates zi and Z2 to zero in (|6l), so that 
in each quark- antiquark operator the two fields are separated by a light-like distance, as they are for 
the usual PDFs. However, the operators describing the first and the second quark are still separated 
by a transverse distance y. The resulting distribution does therefore not involve a twist-four operator 
but rather the product of two operators with twist two. 

For each of the two quark-antiquark operators there are three Dirac matrices Fa that give a leading 
contribution to the cross section, namely Fq = ^7^, F^q — -^7^75 and F^^^ = ^icr"'''"75 with j — 1,2. 
These matrices respectively project on unpolarized, longitudinally polarized and transversely polarized 
quarks and are well-known from the definition of single-parton densities. Note that for multiparton 
distributions one can have spin effects even in an unpolarized proton because the polarizations of 
different partons can be correlated among themselves. If such correlations are large, they leave an 
imprint in the cross section formulae ([5]) and (O , where they can change both the size of the cross section 
and the distribution of particles in the final state. In particular, one finds that a correlation between 
the transverse polarization of two quarks or antiquarks leads to an angular correlation between the 
leptonic decay planes of the two produced gauge bosons. Whether and where spin correlations between 
different partons in the proton are important is therefore an important question for phenomenology at 
the LHC. It is also of great interest from the point of view of hadron structure, and one may hope for 
a fruitful interplay between these two fields in the future. 

Not only the spin but also the color of two partons can be correlated. In Eqs. ([5]) to ([7|) we have 
glossed over this degree of freedom, and it turns out that two-quark distributions can have two color 
structures, both of which contribute to the cross section. The two quark lines with momentum fraction 
xi in Fig. [3^ can couple either to a color singlet or to a color octet. A corresponding color coupling of 
the two quark lines with momentum fraction X2 then follows from color conservation. It is the color 
singlet combination that has the interpretation of a Wigner distribution discussed above, whereas the 
color octet combination describes an interference between different color states of a parton in the 
scattering amplitude and its complex conjugate. The operator defining the color octet combination 
reads (92 ^a2^°92) (9i-fai^°9i), where i° denotes \ times the Gell-Mann matrices. The subscripts 1 
and 2 of the quark fields indicate that after Fourier transformation the fields are associated with 
momentum fraction xi and 2:27 respectively. Although the possibility of spin and color correlations in 
multiparton distributions has been pointed out long ago |14| . the lack of a method to estimate their 
size has prevented them from being included in phenomenological analyses up to now. 

In Sect. [T] we mentioned the resummation of Sudakov logarithms in single hard-scattering processes 
with small transverse momenta. Such logarithms also appear in the double scattering process depicted 
in Fig. [3^. They can be computed by generalizing the method of [H; [l^ from single to double hard 
scattering. An important result is that the leading double logarithms in qf/Q^ are given by the 
product of the corresponding factors for each separate hard-scattering process. Beyond this accuracy, 
the structure becomes more involved since soft gluons can be exchanged between all parton lines in 
Fig. [5^1.. Explicit calculation shows that these effects also lead to a mixing between the color singlet 
and color octet distributions mentioned in the previous paragraph. 



4 Povirer counting: when are multiple interactions suppressed and when are they not? 

When deriving the cross section formula ([S]) we have taken the leading approximation in the small 
parameter A/Q, where the large scale is provided hy ^ q\ ^ whereas the small scale A is 
given by the transverse boson momenta |qi|, |q2| or by the scale of nonperturbative interactions, 
whichever is larger. One readily finds that, up to logarithmic terms, the overall scaling behavior of the 
double-scattering cross section ([S]) is 

Y[l^^dxax^<Pcl, ^ Q^A^ ■ 

One finds the same behavior for the case where the two bosons are produced in a single hard scattering, 
as illustrated in Fig. [3}d. In the fully differential cross section, multiple hard interactions are thus not 
power suppressed with single hard scattering and can hence be of significant size. 

The situation changes if one integrates over the transverse momenta qi and q2. In the double 
scattering mechanism both |qi| and |q2| result from transverse parton momenta and are hence limited 
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to size A. By contrast, with a single hard scattering the sum |qi + q2| is of order A but the individual 
transverse momenta can be as large as kinematically allowed, i.e. they can be of order Q. Single hard 
scattering can thus populate a larger phase space, and one has 

nLi dxi dxi Q'^ ' nLi dxi dxi 

The double scattering contribution is now power suppressed in the overall cross section. This is in- 
deed necessary for the consistency of the usual factorization formulae, which only describe the single 
scattering contribution. 

Multiple interactions are thus important in observables that are sufficiently exclusive and in kine- 
matics where this mechanism can contribute. The relevant kinematic regions are those where one has 
subsets of final-state particles for which the vector sum of transverse momenta is small compared with 
the hard scale in the process. In our case this is the double Drell-Yan process with small transverse mo- 
menta of the produced vector bosons. Another important example emphasized in ,17] is the production 
of two dijet pairs, each of which is approximately back-to-back. 



5 High transverse momentum 

The predictive power of the theory increases if one considers kinematics where the transverse momenta 
(\i are large compared with the scale of nonperturbative interactions (while still being small compared 
with the very large scale Q). At least some of the transverse parton momenta are then large as well, 
and one can compute the transverse-momentum dependent multiparton distributions in terms of a 
perturbative subprocess at scale |qi| and distributions that depend on fewer variables. 

Two types of graphs are of particular importance in this context. In the double ladder graph 
of Fig. 13^, each of the two partons that will enter the collision process acquires a large transverse 
momentum by radiating a parton into the final state. This is a straightforward generalization of the 
corresponding mechanism for single-parton densities at large k. The two ladders in the graph are 
independent, so that this mechanism can contribute for interparton distances y of hadronic size. An 
important finding is that the color factors for ladder graphs favor the case where the parton pairs with 
equal momentum fraction xi or X2 are coupled to color singlets. It remains to be studied quantitatively 
how strong this preference for the color singlet channel is. 

In the graph of Fig. a single parton splits into two partons with high transverse momenta, both 
of which will then take part in a hard collision process. This mechanism is in particular relevant for 
large r, which translates to a small distance y. One finds that the transverse-momentum integrated 
distribution F{xi^y) behaves like 1/y^ at small y due to these splitting graphs. This poses a problem 
of consistency, since the integral over y in the cross section ([7]) is then linearly divergent in y^ at 
short distances. Likewise, one finds that the y integral in ([5]) diverges logarithmically at small y. The 
underlying problem signaled by these divergences is that the graph shown in Fig. [S] can be interpreted 
in two ways. If one identifies the two boxes in the figure as representing the splitting graphs for a 
quark-antiquark distribution, one has a graph for double hard scattering as just discussed. However, 
without the boxes one has a graph for single hard scattering, namely for the fusion of two gluons into 
two electroweak gauge bosons via a quark loop. For a detailed investigation of this loop graph we 




Fig. 4 Graphs where partons with high transverse momentum are generated by parton radiation (a) or by 
the splitting of a single parton into two (b). 
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Fig. 5 Graph for the production of two gauge bosons, where in both colliding protons an initial gluon splits 
into a quark-antiquark pair. The boxes indicate the splitting graph of Fig.[4)D for a quark-antiquark distribution. 

refer to [l^ . A proper theoretical formulation will need to make sure that there is no double counting 
problem for this graph and remove the divergences in the y integration of the double scattering cross 
sections ([5]) and ([7]) . To achieve this in a consistent and practicable way is an outstanding problem. 

6 Approximation by single-parton distributions 

To develop a phenomenology of multiple interactions, one needs a simple ansatz for multiparton distri- 
butions as a starting point. For those distributions that admit an interpretation as Wigner distributions 
it is natural to approximate them by the product of single-parton densities, neglecting any possible 
correlations between the different partons. 

To formalize such an ansatz we insert a complete set of intermediate states J^x between 
the two operators {q2 ra^qT) and {qi Fa^qi) in the definition ([B]) of the two-parton distribution. If we 
assume that this sum is dominated by single-proton states, then we obtain an approximation in terms 
of single-parton distributions. For transverse-momentum integrated distributions, this approximation 
reads 

Fa^,a2[x^,y) J (fh fa^{x2 ,h) fa,{xi ,h + y) , (10) 

where fa{x,h) is the impact parameter distribution of parton a discussed in Sect. [1] This approxi- 
mation, visualized in Fig. [51 is the starting point of most estimates for multiparton scattering in the 
literature, although there is a number of attempts to go beyond it [T9I: [20I: [2l|. 

The corresponding approximation for transverse-momentum dependent distributions has a simpler 
structure if we use the variable r rather than y. We then have 

Fa^^a2{xt,'ki,r) « /Q2(a;2,k2 - ^X2r; ~r) fai{xi,ki - ^xir;r) . (11) 

Here /a(a;,k,r) is a transverse-momentum dependent generalized parton distribution, defined by the 
Fourier transform of a matrix element as in ([T]), where the annihilated quark carries momentum frac- 
tion X and transverse momentum k — ^r and the created quark carries momentum fraction x and 
transverse momentum k + ^r. 




Fig. 6 Visualization of the approximation (|10p of a two-parton distribution in terms of impact parameter 
densities for single partons. In the figure we use that both types of distributions can be represented as squares 
of light-cone wave functions of the proton. 



9 



This strategy of approximation can also be adapted to the color octet distributions described 
in Sect. [31 which cannot be interpreted as Wigner distributions. To this end, we perforin a Fierz 
transformation in color and in spin space, 

(92 ra^t"-q2) {qiFa^t^qi) = -\ {q2 ra^q2) (gi-Taigi) + 5 X! Caia2,&ifc2 (?! (72) (^2 91 ) , (12) 

where the coefficients Cq^ 02,6162 S''^^ straightforward to calculate. The first term on the r.h.s. is a color 
singlet distribution and can be approximated as described above. In the second term, the quark and 
antiquark fields coupled to color singlets are now associated with different longitudinal momentum 
fractions. Repeating the steps that lead to ([TT]) one finds that this term is represented by transverse- 
momentum dependent GPDs in which not only the transverse but also the longitudinal parton momenta 
differ between the created and the annihilated quark. This just means that the skewness variable ^ 
mentioned in Sect. [1] is nonzero in this case. 

It should be emphasized that the above approximations are ad hoc in the sense that we cannot 
quantify to which extent the sum over all intermediate states X in our derivation is actually dominated 
by single-proton states. However these approximations may serve as a starting point, and they can 
give valuable information about aspects that are poorly known otherwise, such as the distribution of 
partons in transverse space and the possible size of color interference effects. 

7 Possible strategies for modeling multiparton distributions 

Our current knowledge of multiparton distributions is very limited, and it should be interesting to 
investigate them in approaches that have successfully been applied to single-parton densities. The 
comparison of predictions obtained in such approaches with upcoming data on multiple interactions 
at LHC could lead to significant progress, with benefits for both sides. 

One possible avenue is to take Mellin moments of two-parton distributions and to study the resulting 
matrix elements in lattice QCD as described in [l3l |. This will naturally give limited information about 
the dependence of the distributions on their momentum fractions Xi, but it has the potential to quantify 
parton correlations in a genuinely nonperturbative framework. 

Constituent quark models, chiral quark models, the Dyson- Schwinger approach and similar meth- 
ods should be suitable to evaluate not only the distributions of single partons but also two-parton 
distributions. As illustrated by simple examples in [Tsj . one may expect two-parton correlations of 
significant size in these approaches. They can also be used to investigate the degree of validity of the 
approximations (ITU)) and (dJ). 

The two avenues just discussed will mostly give information about partons with momentum frac- 
tions Xi above, say, 0.1. In multiple interactions at LHC, typical values of Xi are smaller by at least an 
order of magnitude. One way to investigate this region would be to take the results from quark models 
as an input for scale evolution, which will generate partons with lower Xi by perturbative splitting as 
the scale is increased. Obvious questions to study in this context are to which extent scale evolution 
may wash out two-parton correlations present at low scales, and whether evolution to higher scales 
might improve the quality of the approximation in (|10[) . These would be important lessons, even if one 
has to bear in mind that perturbative evolution alone cannot account for sea quarks and gluons at 
small x, as already mentioned in the beginning of Sect. [TJ An outstanding challenge for the study of 
hadron structure is to find and develop nonperturbative approaches suitable to describe these degrees 
of freedom in a quantitative way. 
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